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Abstract 

Solving the Maximum a Posteriori on Markov Random 
Field, MRF-MAP, is a prevailing method in recent inter- 
active image segmentation tools. Although mathematically 
explicit in its computational targets, and impressive for 
the segmentation quality, MRF-MAP is hard to accomplish 
without the interactive information from users. So it is 
rarely adopted in the automatic style up to today. In this 
paper, we present an automatic image segmentation algo- 
rithm, NegCut, based on the approximation to MRF-MAP. 
First we prove MRF-MAP is NP-hard when the probabilis- 
tic models are unknown, and then present an approximation 
function in the form of minimum cuts on graphs with nega- 
tive weights. Finally, the binary segmentation is taken from 
the largest eigenvector of the target matrix, with a tuned 
version of the Lanczos eigensolver. It is shown competitive 
at the segmentation quality in our experiments. 



1 Introduction 

Image segmentation is an important research field in 
computer vision and computational graphics. In recent 
years, many interactive segmentation methods, such as 
GrabCut([l]), Paint([2]), are designed in the Bayesian 
framework, to say, Maximum a Posteriori (MAP) on 
Markov Random Fields (MAP)([3]). Due to the high quality 
of their segmentation results, especially of accurate bound- 
aries of extracted objects/foregrounds, these methods are 
capable of processing many tedious and time-resuming in- 
teractive tasks in image editing, medical diagnosing, etc. 

In MRF-MAP-based segmentation methods, the sample 
foreground and/or background zones are designated by user 
interactions first, either in the scribble-based or painting- 
based style. The sample data is then adopted to typically 
produce two appropriate probabilistic models for the color 
distributions in the foreground and background zones. Ac- 
companied with some well-designed smoothness terms, the 
data terms appear in the form of the log-hoods of all prob- 
abilities of pixel colors, then compose the final target func- 



tion. The target function is then to be optimized to obtain 
a good segmentation result. Clearly, the user interactions 
are important in the whole task, in other word, it's the first 
power to push the optimization to run. Especially when the 
smoothness term only reflects the noisy tolerance capabil- 
ity, as in GrabCut and Paint, the user-defined color distribu- 
tions are much critical in determining which part the pixels 
belong to. 

Compared with the generally dissatisfactory segmenta- 
tion quality in current automatic segmentation methods, it's 
tempting of considering the good segmentation results in 
MRF-MAP-based interactive tools. So an interesting prob- 
lem arises naturally here: whether or not we can design a 
segmentation algorithm in the MRF-MAP framework, but 
without any user interactions? However there are few con- 
crete attempts except SE-Cut, to the best of our knowledge. 
In order to produce the required color models, the user- 
designated zones in GrabCut are replaced by automatically 
chosen seeds in SE-Cut([4]). SE-Cut chooses its seeds by 
random walks and spectral embedding techniques instead in 
the MRF-MAP framework itself. Clearly, it's not a unified 
processing in only one dominating thread in the theoretical 
sense. In practice, there's an additional but necessary re- 
quirement for all seeds to be large enough lest being elim- 
inated in the advanced processing. However, it's hard to 
derive from the MRF-MAP directly, and it makes SE-Cut 
a little more inconvenient both in practice and theoretical 
views. 

In this paper, we present an automatic segmentation al- 
gorithm, NegCut, thoroughly under the MRF-MAP frame- 
work. It is designed to produce binary segmentation re- 
sults, and likely to figure the dominating objects or out- 
standing scenes out of the images, the same as the inter- 
active segmentation tools. Our main contributions are: (1) 
a mathematically strict analysis on the hardness of MRF- 
MAP without any predefined information (or uninformed 
MRF-MAP as abbr.) (2) a feasible way to approximate the 
MRF-MAP functions, and a concise mapping from the op- 
timization of approximating energy onto minimum cuts of 
graphs with negative weights. (3) a simple implementation 
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with the help of the Lanczos eigensolver([5]), along with 
some experiments and corresponding analysis. 

2 Background 

In a MRF-MAP-based interactive segmentation algo- 
rithm, the task is mapped into the minimization of a par- 
ticular target energy function E, which is the sum of a data 
term Erj and a smoothness term Eg multiplied by a factor 
A: 

E = E D + XE S (1) 

where the data term Ed is the sum of all data penalty on 
all pixels according to their labels, the smoothness term 
Es is the sum of all potential of adjacency interactions be- 
tween all neighboring pixels of different labels. In our sit- 
uation, these labels are just foreground or background (fore 
and back as abbr.) Typically in current algorithms, the data 
penalty of each pixel is the negative log-hood of its proba- 
bility according to the color distribution in the foreground 
or background zone, and the smoothness terms are usually 
defined on all adjacent pixels on 4/8 -connection grids: 

E D =J2-^Pr L(p) (p),E s = J2 S &l) 

P (p,q)£NMp)^L(q) 

(2) 

where L(p) is the label of pixel p in the label configuration 
L, and (p, q) <G N means p is adjacent to q. Although the 
smoothness terms are critical and important also, it's not the 
point in this paper. In general the smoothness term could 
be regarded as a prior in form of a Markov Random Field, 
therefore the problem is actually to pursue an appropriate 
label configuration to reach the lowest energy E, which is 
the Maximum a posterior under the MRF-MAP framework. 
Although it's more exact to regard the prior as a conditional 
random field or discriminating random field, we here still 
take MRF as the notation following the past literatures for 
convenience. 

There are many choices when determining the coeffi- 
cients in the energy function E. For the data term, the or- 
dinary statistical histograms are sufficient for the situations 
with only a few possible colors, such as 16 gray levels, 256 
colors and so on. However, it is not suitable for large color 
spaces such as 24-bit true colors, as the samples were rela- 
tively too sparse to raise a meaningful probabilistic model. 
Then other more complicated models, such as the Gaussian 
mixture model (GMM), fit this situation better considering 
the compromise both of efficiency and accuracy. It should 
be pointed out that, the GMM method is compatible to the 
histogram method, if we take each color value as an exact 
Gaussian component with variance 0. 

For the smoothness term, the coefficients should be de- 
signed to encourage continuity between adjacent pixels, es- 
pecially those of similar colors. A good choice is the ex- 
ponential function of the Euclidean distances of colors, as 



already adopted in GrabCut. Also we could add a posis- 
stive constant to all smoothness coefficients to enhance the 
local continuity, or reduce the factor A down to eliminate 
the differences between different colors. 

After the color probabilistic models and smoothness co- 
efficients are determined, there are several methods to min- 
imize the target energy function, such as graph cut, loopy 
belief propagation (LBP) and iterated conditional mode 
(ICM)([6]). In the case of binary labels, the graph cut 
method is shown to be practically efficient: the target energy 
function could be mapped into cuts on an undirected graph 
where the global minimum solution could be obtained with 
the max-flow algorithms in the polynomial time. Here the 
data penalties are regarded as the weights of edges linked to 
the source/sink node in networks. As to the case of multiple 
labels, however, the hardness changes, to say, it's a k-way 
cut problem which is proven NP-hard. 

The determined coefficients of the data terms, derived 
from the probabilistic color models, are important to con- 
struct the graph to be cut. Undoubtedly the user interac- 
tions are necessary for the current MRF-MAP based in- 
teractive image segmentations, whereas other tools such 
as random walks could present some candidate fore- 
ground/background zones instead. 

3 Hardness of uninformed MRF-MAP 

Remind that in the interactive image segmentations 
based on MRF-MAP, the probabilistic color models are 
computed from sample pixel colors, to be more detailed, 
foreground/background zones designated by users. It means 
that the whole energy function is still computable even if we 
only know the exact labels of all pixels. From this point of 
view, now the energy function E can be written as 

E(L)=Y,-^Gen(L,p) + \ ]T S(p,q) 

P (p,q)£N,L(p)^L(q) 

(3) 

where we aim to obtain an appropriate probabilistic color 
model from all pixel p's with the same label L p in {fore, 
back}, with the help of a histogram- or GMM-generating 
function Gen. Given Gen, all we need to calculate E is just 
the labels of all pixels in the image. In other word, there is a 
corresponding energy value for each label configuration of 
pixels. Clearly there exists a minimum energy E* according 
to a certain label configuration L* among them: 



V 



arg min(£'(L)) 
L 



(4) 



Recall that we usually obtain a good segmentation result 
according to the minimum energy value in the interactive 
segmentations, then whether or not L* is also meaningful 
in the image segmentation problem? Furthermore, could we 
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obtain such a global optimal solution L* in the polynomial 
time? However, it is not addressed before to the best of our 
knowledge. 

Since the histograms are easy to deal with, it's reason- 
able to investigate the hardness of uninformed MRF-MAP 
on small number of colors, while the conclusion still holds 
under the general cases of GMMs. To be more concen- 
trated, we neglect all smoothness term, i.e., A = 0. Sup- 
pose we have n pixels of m colors in the input image, and 
the pixels of color i counts to rij. Then for each label 
configuration L in {fore, back}", we set no,, as the pixel 
amount of color i and label fore, while it's the same for 
ni : i to color i and back. Furthermore, F and B are the 
two sets of fore pixels and back pixels respectively, while 
their pixel amounts are so = no,i + no. 2 + • • • + no,m an d 
s% = Hi i + "1.2 + . . . + ni, m . These notations will also be 
taken in the following sections. Now, the energy function E 
on L becomes 



p£F U pEF 

m 

= (s In s + s\ lnsi) - ^(n ,j lnn ,j + n lti lnn M ) 

i=l 

(5) 

Then consider a function f\ defined on interval [0, c]: 

/i (0) =a In a + (b + c) ln(6 + c) - c In c, 
/i (c) =(a + c) ln(a + c) + b In b — c In c, 
fi(x) =(a + x) ln(a + x) + (b + c — x) ln(6 + c — x) 
— xhix — (c — x) ln(c — x), (0 < x < c). 

(6) 

Obviously /i is continuous on the entire interval, and we 
have its derivative f[(x) = on x = f[{x) > when 
< x < and f[(x) < for ^ < x < c. Hence 
the maximum value of /i on [0, c] appears on x = 
whereas its minimum value is on x = or x = c. Then for 
each k = 1, . . . , m, if sq, s\, and all uqj, ni t i(i ^ k) as 
fixed, we would get the minimum value of E(L) at no,fc= 
or fii fc= 0. Without the loss of generality, we assume there 
exists an m* satisfying that n\ t k= for all k < m* , and 
n>o,k= for all k > m*. So we have 



k<r 



k<r 



k>m* 



iE(L) =( n k )H Yl »*) + (E nfc ) ln ( E n ^ 

m 

E n fc m n k 



k=l 



(7) 
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Figure 1. An 1 x 2k image with m colors. 

that /2 reaches its minimum value at x = |. So finally we 
know that the energy function E reaches its minimum value 
when^,,^, rik = Ylk>n* Uk ^ possible. In the following 
we will reduce the set partition problem to the minimization 
of our energy function E. Since the set partition problem is 
an NPC problem([7]), uninformed MRF-MAP is NP-hard 
then. Remember that, the set partition problem is about 

Given a positive integer set X ={x\, X2, ■ ■ ■ , x m } sum- 
ming up to 2k, whether or not there is a subset X' satisfying 
that the sum of all entries in X' equals k? 

The reduction is rather straightforward. First we con- 
struct an image img with the size 1 x 2k (only one pixel 
high) in the color space {color i, color 2, ■ ■■ , color m } as 

Then in the smoothness terms, we make it sufficiently 
large for pixels of the same colors, but close to zero for 
pixels of different colors. Therefore the blocks of identical 
colors would never be cut, and the extremely small coeffi- 
cients for adjacent pixels of different colors would have no 
effect on the energy function E. When minimizing E, we 
would get 




Now consider another function f2(x) = x\nx + (d — 
x) \n(d — x). Since /^(x) = In ^3^, it's easy to conclude 



on a certain subset X' satisfying 2~2ieX' ^ l ~ Z-^iex r 
k, according to the above analysis. So the set partition 
problem could be determined by minimizing the energy 
function E and verifying whether the minimum is equal to 
(2k In k — YhLi x i m x i)- So the reduction is finished. 

Here we have established the link between the set parti- 
tion problem and the uninformed MRF-MAP in image seg- 
mentations. Recall that, the hardness of normalized cut is 
also linked to the same NPC problem([8]), and the reduction 
there implicitly ensures that its binary segmentation results 
are somewhat fair in sizes. Clearly this favorable property 
also holds for the potential image segmentation algorithms 
based on uninformed MRF-MAP, to say, no much isolated 
small-sized pieces were produced. 

4 Approximating uninformed MRF-MAP 

Considering the NP-hardness when minimizing the en- 
ergy function E, it is reasonable to pursue its approximating 
solution instead. However, E is not a polynomial with re- 
spect of the label configuration L, even not a closed form at 
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all. So it is feasible to find its replacement which is solvable 
and sufficiently close to E. 

First let's consider a function/3 (x) = x\nx + {1 — 
x) ln(l — x). We have fz(x) = ^xy — A(x) where y 
= (1 - x) and 



A(x) = xy(^{x 2 + y 2 ) + -{x 3 + y 3 ) + 



), (9) 



by performing two Taylor expansions. The mean of A(x) 
on [0, 1] is f Q A(x)dx = and we could basically taken 
it as the value of A(x) on the entire interval, given that 
the mean square error L (A(x) — j^) 2 dx ~ 3 x 1CP 4 is 
considerable small. Hence it is totally acceptable to ap- 
proximate /3 with — x) — yg), as shown in Fig. 2. 
Now return to our energy function E. Since the uninformed 




Figure 2. Curves of function /3 and its approximation on [0, 1]. 

MRF-MAP on small number of colors is relatively simple, 
we set the case as our start. 

4.1 For small color spaces 

Histograms are generally adopted in the cases of small 
color amounts. Based the approximation of f$, we can ap- 
proximate the energy function E with 

(^a (n - so) " ^) " g "<m("' ~ ~ 1^) 



A E s (p>q) 

(,p,q)eN,L{p)^L(q) 



(10) 



Since m, n and all m's are all fixed in the input image, it is 
actually to minimize 

—5 m —5 

— S (n - 3o) - E n 0,i( n i - n 0,i)) 

1=1 (11) 
+ A S M 

(p,q)€N,L(p)^L(q) 

when pursuing the minimum of E. 



4.2 For large color spaces 

It is a little more complicated for the cases of large color 
amounts, such as 2 24 colors in the 3-bytes RGB space. 
Since the amount of potential colors is usually much larger 
than that of the image pixels, there is hardly remarkable 
number of samples in each bin of the color histogram. 
GrabCut takes GMM as a good replacement to histograms 
in the color image segmentation, but it is not fluent and easy 
to find the approximation of E from GMMs. Instead we 
adopt a trivial but also accurate scheme in the continuous 
RGB color space: first, choose a fidelity parameter a 2 em- 
pirically or by calculations on the entire image as in [1], 
then set the statistical contribution of sample color on each 
pixel p to be a normal distribution with mean p and variance 



Pr p {k) 



1 



V2~7 



(12) 



where Dist(p, k) means the squared Euclidean distance be- 
tween the colors of pixel p and k. Then construct the prob- 
abilistic color models on all foreground pixels and back- 
ground pixels as 

Pr F (p) = — E Pr q (p),Pr B {p) = — V] Pr q ( P ) 
rig * — ' m * — ' 



q£F 



qeB 



(13) 

And the energy function E becomes into the smoothness 
term plus an integral on the entire continuous RGB space 

E(L) = J (-^ • Pr F (p) In Pr F (pj) dp 

+ J (-^ ' Pr B {p)\nPr B (p)j dp (14) 

(p,q)eN,L(p)^L(q) 

Since^ ? Pi q (p) is a fixed value in the input image, it is 
actually to minimize 



W__5_ 5 f Pr p (fc) ■ Pr g (fc) 
2ni 2n J ^Pr^fc) 

(p,q)&N,L{p)jtL(q) 



dk 



(15) 



when pursuing the minimum of E. Again we can make it 
simpler by enlarging the factor A n times 



+ A J2 s &>ti 

(p,q)£N,L(p)^L(q) 



(16) 
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Figure 3. Leftmost column: original gray images, other three columns: segmentation results with A = 1, 5, 10 in turn. (16gr ay scales) 



5 NegCut based on Min-Cut 

5.1 NegCut on small color spaces 

Take the simple case of histograms on small number of 
colors as our start. Then our target is to obtain a segmenta- 
tion result by minimizing (11). Now construct an undirected 
graph G of n nodes corresponding to the pixels, and set the 
edge weight w(p, q) as the sum of the following three terms 

-5 
2n 

■J*- , p and q have the same color i, 
0,ofherwise. 

S(p, q),pis adjacent to q, 
0,otherwise. 



w\{p,q) 
wa(p, q) 
w 3 {p, q) 



(17) 



It is easy to prove that the approximating energy value on 
a label configuration L is just equal to the capacity of the 
corresponding cut C= {F, B} on G. Therefore it is to find 
the minimum cut on G minimizing (11). Although there 
exists several 0(N 3 ) algorithm solving the minimum cuts 
on graphs with non-negative weights, it is not suitable here. 
In fact, the negative weighted edges make the Min-Cut to be 
a Max-Cut on a auxiliary graph with non-negative weights. 
Since Max-Cut has already be proven to be a NP-complete 
problem, it can be concluded that there is no polynomial 
algorithm for us to minimize (11). 

However, we might generalize this discrete optimization 
task into the continuous real space R n . First we set an in- 
dicator vector D = [di, d%, ■ ■ ■ , d n ] T e {+1, —1}" denot- 
ing the ith pixel of fore by +1 or back by —1, then estab- 
lish a matrix W = [w(p, q)] whose entries sum up to Sw- 
Clearly the cut value is just equal to h (Sw — D T WD). 
Subsequently, we generalize the labels into the continuous 
interval [— 1, 1] instead of {+1, —1}, and the original opti- 
mization task on (1 1) becomes 



According to the Lagrange factor method([6]), the solution 
to (18) is the eigenvector D corresponding to the largest 
eigenvalue. Then we might straightforwardly get the re- 
quired binary indicator from D by setting the ith entry as 
— 1 if di < 0, or +1 if di > 0. Here the Lanczos algorithm 
is adopted to finish the time-consuming operations of the 
largest eigenvectors for large matrices. In fact, the Lanczos 
algorithm is well known as the fastest method to solve the 
largest eigenvectors in O(N) time mostly for sparse matri- 
ces, whereas our graph G is completely connected and the 
corresponding weight matrix W is full. Fortunately, W is 
only called to perform a matrix- vector multiplication there, 
and the major calculation is identical for pixels of the same 
color. So we might improve it to sharply reduce the time 
consumption down to linear, as described in algorithm 1 . 

Now we have established the NegCut algorithm on the 
small color spaces as below: 

Algorithm 1 NegCut on small color spaces 

1: On each adjacent pair of pixels (p, q), compute the 
smoothness term S(p, q)\ 

2: Solve the largest eigenvector D = [d± , c?2, • • • , d n ] of 
the matrix W = [w(p, q)] with the Lanczos algorithm, 
but perform the matrix-vector multiplication of W and 



some vector R = [n , ra, ■ • 
• set ip = 0. Then for k ■■ 



, r„ as 



1 ton, cp <- ip - ^ • r k ; 



m&xD T WD,s.t IIDII 



(18) 



• set 0i = for all is. Then for k = 1 to n, set 
d L 0i + ■ rfc according to its color i\ 

• for k = 1 to n, calculate [i k = J2j:( 3 ,k)eN ^0' 
Then according to its color i, and output <p + di + 
( In ~ 2n" ) ' Vk as tne ^ entr y °f vector W ■ R. 

3: let the label of pixel k be back if dk < 0, or fore if 
d k > 0. 
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Apparently we could finish step 1 and 3 in 0(n) time. 
From step 2. a to 2.d, it requires 0(n) time to calculate the 
product vector involving W, so we need totally 0(n) time 
for the slightly changed Lanczos algorithm in step 2. Finally 
the entire time complexity of NegCut for small color spaces 
remains 0(n). 

5.2 NegCut on large color spaces 

The case of large color spaces is a little more compli- 
cated. Recall that in (16) we need to calculate the integral 
on the whole continuous color space, and it is unrealistic in 
most applications. A feasible substitution is to perform the 
calculation on the typical samples according to the prob- 
abilistic distribution of (13). Fortunately the pixels of the 
input image just meet this requirement well. Hence the in- 



Pr p (fc)-Pr,(fc) 



however 



tegral becomes into the sum as J2k y\ p r (fc) 
the computational load is still heavy as we have to scan each 
pixel and record its contributions to n Pr 9 's. While the en- 
tire complexity is up to 0(n 2 ), which is totally unaccept- 
able even on ordinary-sized images. In order to improve it, 
we cluster all colors into a limited number of classes first, 
so that the summations would be identical for the pixels in 
the same color class. Now, all we need is just to scan each 
pixel and record its contributions to these color classes, and 
the time consumption is reduced to 0(n) given the amount 
of color classes is limited. 

As in the small color spaces, we construct an undirected 
graph G of n nodes corresponding to the pixels, and set the 
edge weight w(p, q) as the sum of following three terms, 



w 2 {p,q) = 2 ' 

w 3 (p,q) = 



Pr p (fc) • Pv q (k) 



S(p, q),pis adjacent to q, 
(Xotherwise. 



(19) 



Now we have established the NegCut algorithm for the large 
color spaces: 

There are some light differences between the two algo- 
rithms in step 1 and 2. However given the amount of color 
classes is limited and fixed, it requires also O(n) time to fin- 
ish the multiplication involving W. Finally the entire time 
complexity of NegCut for large color spaces is also O(n). 

6 Experiments 

Since NegCut is different on the small color spaces and 
large color spaces, our experiments are finished on both 
gray scale images and color images. We trivially set all 
smoothness terms to be 1 over entire images, and let A vary 
from 1 to 10. All color test images are chosen from the seg- 
mentation datasets from Berkeley and MSR at Cambridge, 



Algorithm 2 NegCut on large color spaces 

1: Cluster all pixels into m color classes with the mean 
color ci, . . . , c m . On each pixel k, calculate its prob- 
abilistic densities in all m Gaussian Distributions, and 
add its contribution Pl ^ fc p i Pl ^ fc ^ to w 2 (p, q) for every 
two color classes of pixel p and q; 

2: Solve the largest eigenvector D = [di, cfe, • " > d n } T of 
the matrix W = [w(p, q)] with the Lanczos algorithm, 
but perform the matrix-vector multiplication of W and 
some vector R = [n, r 2 , • • ■ , r n ] T as 

• set ip = 0. Then for k = 1 to n, tp <— tp — ■ r&; 

• set 8i = for all i's. Then for k = 1 to n, i = 1 

to m, set 6i <— 6i + Wz(j, k) ■ r k according to its 
color class j; 

• for k = 1 to n, calculate /j, k = J2j : (j,k)&N S U> k )- 
Then according to its color class i, output ip + L + 



(J-k 



2n 



r k as the fcth entry of vector W ■ R. 



3: let the label of pixel k be back if dk < 0, or fore if 
d k > 0. 



and all gray scale images are derived from them. Here the 
amounts of color classes in all color images are set to be 16, 
and all gray-scaled images also have 16 gray levels within 
them. Based on the block-based Lanczos eigensolver([9]), 
we established two versions for NegCut in Matlab codes. 
At last, all test images are resized to be about 256 x 256 to 
reduce the time-consuming calculations in the experiments. 

There are three groups of results in Fig. 3, 4 and 5: gray- 
scaled on different A's, segmentation results of color im- 
ages on different A's, and different color images on fixed 
A = 1. Our target images are intendly chosen to be of close 
gray scales, splendidly colors and delicate local changes, so 
that to verify the performance of NegCut on different chal- 
lenges. The small color space version of NegCut works on 
gray-scaled images, while those color images are processed 
with the more complicated algorithm. 

In general, we obtained basically acceptable segmenta- 
tion results in all experiments, especially that most domi- 
nating objects are figured out of these images. There are 
more isolated, but vivid pieces in the segmentation results 
when A = 1, both on gray-scaled images and color images. 
And it's much better for color images because the connec- 
tions between different colors are involved in the large color 
space version of NegCut. When A = 1, the segmentation 
boundaries are more likely located on the desired edges of 
the objects in these images. However, it also brings too 
much emphasis on these discontinuous line segments, and 
results in much more isolated pieces in the segmentations. 
When A varies from 1 to 5, then to 10, it is shown the seg- 
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Figure 4. Leftmost column: original images, other three columns: segmentation results with A = 1, 5, 10 in turn, (m = 16) 



mentation boundaries become smoother and smoother, but 
on the cost of the lose of elaborate details. When A = 10, 
the two segmentation zones are more fair in sizes. The rea- 
son for that is, the large portion of the smoothness terms 
in the energy function weakens the effect of the data penal- 
ties. NegCut tends to cut two pieces sized of one half of the 
entire image to reach the minimum energies. 

7 Conclusions 

Despite the NP-hardness of uninformed MRF-MAP, 
NegCut is an initial attempt in developing automatic image 
segmentation tools under the MRF-MAP framework. Given 
that the Lanczos algorithm is specially designed for solv- 
ing the extreme eigenvectors, NegCut are much easier to 
develop in principle than those require inner eigenvectors, 
such as the normalized cut. Though NegCut only produce 
binary segmenation results, it can be recurrently called to 
refine the past segmentationzones results until the expected 



one appears. 

Recall that in the interactive segmentation tools, their 
target is to obtain the minimum of particular energy func- 
tions with user interactions. Since NegCut approximates 
the global minimum energy value, it is meaningful in the 
analysis on the underlying mathematical trends when these 
interactive tools performing calculations. 

Different from the ordinary minimum-cut-based seg- 
mentation or clustering algorithms, the negative weights 
ensure that NegCut likely obtain fair segmentation re- 
sults instead of isolated extremely small pieces. However, 
the double-fold NP-hardness encountered in solving MRF- 
MAP reminds us, it is necessary to make a attempt on more 
efficient calculation methods. 
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